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Abstract 

It has been observed that a quantum mechanical theory need not to be Hermitian to 
have a real spectrum. In this paper we obtain the eigenvalues of a Dirac charged particle in 
a complex static and spherically symmetric potential. Furthermore, we study the Complex 
Morse and complex Coulomb potentials. 
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1 Introduction 



The first interest study in the non-Hermitian quantum theory date back to an old paper by 
Caliceti et al [1]. In [1] the imaginary cubic oscillator problem in the context of perturbation 
theory has been studied. The energy spectrum of that model is real and discrete. It shows that 
one may construct many new Hamiltonians that have real spectrum, although, their Hamiltoni- 
ans are not Hcrmitian. The key idea of the new formalism (non-Hcrmitian quantum theory) lies 
in the empirical observation that the existence of the real spectrum need not to necessarily be 
attributed to the Hermiticity of the Hamiltonian. In non-Hcrmitian Hamiltonian, the current 
Hcrmiticity assumption W = H is replaced by the PT-symmetry condition as = PTHPT 
[2, 17], where P denotes the parity operator; (Pip)(x) — ij)(—x) and T the time reversal op- 
erator; (Tip)(x) = ip*{x) or generally by Pseudo-Hcrmitian condition as W = t~ x Ht, where 
r is an invertible Hermitian operator [17]. Such non-Hermitian formalism, for the context of 
Schrodingcr Hamiltonian, has been studied for many different subjects with several techniques 
[1-20]. Also, some explicit studies of the Hermitian and non-Hcrmitian Hamiltonians have 
performed in the context of Dirac Hamiltonian. For example, the solution of ordinary (Hcr- 
mitian) Dirac equation for Coulomb potential including its relativistic bound state spectrum 
and wave function was investigated in [22,23]. Also by adding off-diagonal real linear radial 
term to the ordinary Dirac operator, the relativistic Dirac equation with oscillator potential 
has been introduced [24, 25] moreover the energy spectrum of corresponding eigenfunctions 
have been obtaind. The ordinary (Hcrmitian) Dirac equation for a charged particle in static 
electromagnetic field, is studied for Morse potential [26]. 

In this paper, we consider the non-Hcrmitian Dirac Hamiltonian for complex Morse and 
complex Coulomb potential. We consider a charged particle in static and spherically symmetric 
four component complex potential. By applying a unitary transformation to Dirac equation, 
we obtain the second order Schrodingcr like equation, therefore comparison with well-known 
non- relativistic problems is transparent. Using, correspondence between parameters of the two 
problems (the Schrodinger equation and the Schrodinger like equation which is obtained after 
applying the unitary transformation on Dirac equation for a potential ) we can obtain the 
bound states spectrum and wave function. 

The structure of this article is as follows. In sec. 2, we study the non-Hcrmitian version of 
Dirac equation for a charged particle with static and spherically symmetric potential, then by 
applying a unitary transformation we obtain the proper gauge fixing condition and Schrodinger 
like differential equation. In sec. 3, we discus the Dirac equation for complex Morse potential 
then we obtain the real energy spectrum and corresponding eigenfunctions. In sec. 4, we consider 
the Dirac equation for a complex Coulomb potential and we obtain the real eigenvalue of it. 

2 Preliminaries 

The Hamiltonian of a Dirac particle for a complex electromagnetic field is (c = % = 1) 

H = a.{p-eA)+f3m + eV, (1) 

where the Dirac matrices a, (3 have their usual meaning, and setting Aq equal to V. In (1), A 
and V are the vector and scalar complex field respectively, where A* ^ A and V* ^ V . Then 
the Dirac Hamiltonian (1) is not Hermitian. It is well known that the local gauge symmetry in 
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quantum electrodynamic implies an invariance under the transformation as: 

(V,A)^(V',A') = (V + ^,A + ^A). 



(2) 



Here A(£, ~?) is a complex scalar field. Suppose that the charge distribution is static with 
spherical symmetry, so the gauge invariance implies that V — V and A 1 — fA(r) , where r is 
the radial unit vector[26]. One can denoted the correspondence wave function of (1) as: 



* = 

In this case one can obtain 

(m + eV - E r )$ = i 

(eV — m — E r )x = i 



(3) 



a.^ — e(a.f)A(r 
a.^ + e(a.f)A(r 



X, 



(4) 



Here it's are the three Pauli spin matrices, E r is the relativistic energy eigenvalue, then we 
replaced iea.A(—iea.A) in first(second) equation of (4) instead of ea.A, respectively. Note 
that, because of the spherical symmetry of the complex field, the angular- momentum operator 
J and the parity operator, P, commute with the Hamiltonian and the two spinors $ and \ 
have opposite parity also. So the correspondence wave functions are denoted by 



It is seen that 



X = /(^t!J- v (^). 

. Ik 
(a.V)ig(r)n K ^('d,(p) = ia r Q, K ^{d r H \- -)g(r), 



)f(r), 



where k is the spin orbit coupling operator which defined as: 

k = a.L + HI. 

and we have used from 



in which 



kn TK .^(i3, if) = ±Kfi,Q TKt n('&,<f), 

-(/ + !) = -(j+i) for j = l + \ 
l = U + b for j = l-l 



(5) 

(6) 
(7) 

(8) 
(9) 
(10) 



Therefore by defining u\ — g{r)/r , U2 = f(r)/r, we obtain the following two component radial 
Dirac equation [28] 

(m + eV - E r ) Ul (r) = {8 r - - - eA{r))u 2 {r), 



(eV — m — E r )u2{r) 



(d r + - +eA(r))m. 



(11) 
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Note that, A(r) is a gauge field, which has a symmetry such as (2), therefore, it must be fixed. 
It is seen that fixing this gauge degree of freedom by .~A* = Q£ — is not a suitable choice. 
Remark that in this paper, instead of solving Dirac equation we want to solve the 2 th order 
differential equation, which is obtained by eliminating one component of equation (11). Note 
that, this result is not Schrodinger like. One can obtain the proper gauge fixing by applying 
the global unitary transformation on two components u\ and u 2 as: 

Mi = a<f> u - b(f> 1 , 

u 2 = b^ u + b^ 1 , (12) 

where a, b £ 3?, a 2 + b 2 = 1 and (j> u ,<j> 1 are the upper and lower component of spinor. This 
unitary transformation create two results. Firstly, it makes a gauge fixing condition such as: 

eV(r) = S(eA(r) + -), (13) 
r 

where, S= 2ab. Secondly, the final 2 th order differential equation is Schrodinger like equation. 
So by applying (12) on (11) and use of (13), we have: 

(2eV + Cm-E r )<i> u +( < ^-Sm-4-)^ 1 = 0, (14) 
\ o ar J 

^-Sm+4-)r + (-Cm-E r )^ 1 = 0, (15) 
S ar J 

where C — a 2 — b 2 . However, by eliminating (f> 1 in (14) and (15), one can obtain the Schrodinger 
like differential equation for radial upper component, (f> u , as: 

-^ + V eff r + (m 2 -E 2 )r=0, (16) 

where 

V eff ^(fv) 2 + (2eE r V-§^). (17) 

3 The complex Morse potential 

The complex Morse potential in Schrodinger equation, which holds discrete spectrum is given 
by [18] 

V CM (x) = (B R + iBtfe- 2 * - (B R + + 2D)e~ x , (18) 

the corresponding Schrodinger equation is as: 

-^ + V CM (x)yj~E^ = 0. (19) 

This Schrodinger equation is exactly solvable. It is well known that the energy eigenvalues of 
the Schrodinger equation for a complex Morse potential (18) is a function of (D + 5 ), which 
explains that the spectrum of that is real [19], which is obtained as: 

E n = -(n-D) 2 , (20) 



4 



where n is an integer number between and D. So we assume that, the complex potential V(r) 
in (17) is: 

V(r) = -(( + i V )e- r . (21) 

Here £,77 £ Jt . By using equation (17), one finds the effective complex Morse potential, V^jy , 
as: 

Vfff = C-§) 2 (C + «?) 2 e - 2r - ^(C + W - 2eE r )e- r . (22) 

By comparing equations (16) and (22) with (19) and (18), respectively, we can obtain a corre- 
spondence between those parameters as: 



c 



Ce 

~S 

Ce 

-2eE r 
m 2 -El 



Br, 

Bi, 
2D, 
-E. 



(23) 



Therefore, from (20) and (23), it is clearly shown that the real spectrum of Dirac particle with 
complex Morse potential is: 



Er, n = { - en + V^l + e 2 )™ 2 -n 2 j 



(24) 



where n = 0, 1, 2, ....,n max < m\f\ + e 2 . Lastly, we can obtain the eigenfunction of (16), 4> u (r), 
with complex effective potential (22) as: 



C( 2 )=^/"e-5^«(4 
where, Af n is the normalization constant and 

z = 2(^)(C + ^)e-, 



(25) 

(26) 
(27) 



One can show that for the case n < —eE r ^ n , the </>^(+oo) = 0, and the orthogonality condition 
is satisfied [19]. Furthermore, from (15), one can obtain the lower spinor component as: 



mC + E r . n 



(Sm + u n + n)Li^(z) + (2u n + n)L 2 ^ 1} (z) 



(28) 



(29) 



where (15), (25), (26), (27) and recursion properties of the Laguerre polynomials have been 
used. 
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3.1 The complex Coulomb potential 

We assume that the complex Coulomb potential is: 

iZa 
eV = , 



(30) 



where a is fine structure constant and Z is atomic number. This equation represents the 
interaction of a point charge, -e, and an imaginary point charge iZe. It is obvious that one 
has the freedom to use any unitary transformation. So for two coupled first order differential 
equation (11), with complex Coulomb potential (30), is obtained the gauge fixing condition and 
Schrodinger like differential equation, by applying the another global unitary transformation, 
instead of (12). In this case we apply the global unitary transformation such e^ 6 ' <Ta= , which is 
rewritten as: 

a ib 
ib a 



U = 



(31) 



where a, b £ 3? and a 2 + b 2 = 1. By applying (31 ) to <f> u and <j> 1 and institute it in (11), we have 

~i{S 2 -C 2 ) 



(m - E r C)<j) u + 
i{S 2 - C 2 ) 



V - iSE r + d r 



-eV - iSE r - d r 



(to + E r C)(j>- 



' = 0, 

= o, 



where, S — 2ab, C = a 2 — b and we have used from a gauge fixing condition as: 

eV-§(eA+J). 



(32) 



(33) 



However, we eliminate the 4> l component in (32), and obtain the radial differential equation for 

4> u as: 

+ Vef C f(rW - m 2 W(r) = 0, 



where 



in which 



KS?(r) = 



7(7 + 1) 2iaZ(S 2 - C 2 )E r 



(34) 
(35) 



(ff 2 - C 2 )aZ 
S ' 
Vn 2 + a 2 Z 2 , 



(36) 



(j + -) 2 + a 2 Z 2 , 



is the relativistic angular momentum. Since the wave function has to be normalized, we have 
chosen the positive sign for 7. Furthermore, it is obviously seen that these solutions, in fact, 



6 



exist for all values of Z 2 . The non Hermitian and PT-symmetric radial Schrodingcr-Coulomb 
differential equation is [10] 



l_ + 1(1 + 1) + iA_ E 
dr 2 r 2 r 



iP(r) = 0, 



(37) 



and its non-relativistic energy spectrum is: 



E n — 



An 2 ' 

n r + 7 + 1, 
0,1,2,... 



(38) 



So comparing equations (34) and (35) with (37), gives the following correspondences between 
the parameters of two problems for the regular solution of (f) u as : 



7 = I, 
2aZ(S 2 - C 2 )E r = A, 
(E 2 -m 2 ) = E. 



(39) 



Using (38) and (39) we obtain the relativistic energy spectrum for complex Coulomb potential 
as: 

a 2 Z 2 (S 2 -C 2 ) 2 ' 



E r .n = m 



1 - 



(n r + 7+ l) 2 

where n r = n — j + \ is the radial quantum number, n the principle quantum number, 7 
is the angular momentum quantum number, and consequently, 



E r 



m 



2\2 



Z 2 a 2 (S 2 -C 2 ) 



[n-3-\ + yj(.3 + \) 2 + Z 2 a 2 Y 



n = 1,2,3,... 
. 13 5 

3 ~ 2'2'2'"" 

and the upper component of the radial spinor eigenfunction is: 

4> u n {r)=Ny^e-iL%+\r). 



(40) 



j + h 



(41) 



(42) 



(43) 



Here r = x — i9, where x is real(x € [0, 00)), 9 > and W„ is normalization constant [10]. 
It shows that the integration path has been shifted down from the positive. From (35), it is 
understood that the potential is zero in two cases. The first one is Z = and S ^ C, and the 
second one is Z ^ 0, S = C. It is easily seen that for vanishing potential the energy eigenvalue 
is m. Note, for a = cos(^) and b = sin(|), S = 2ab = sin(^) = a 2 — b 2 = cos(^). Hence for 

2 For the ordinary Coulomb potential 7 = Vfc 2 — Z' 2 a 2 and, therefore, one can conclude from 7 that for states 
with k 2 = 1 only solution up to Z ~ 137 can be constructed. For (Za 2 ) > k 2 , in general, the real part of the 



wave function shows an oscillatory behavior and for the states that n 
[28]. 



k, the energy spectrum is imaginary 
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the case which the unitary transformation is e^ %(7m the energy eigenvalue is m. Therefore, we 
assume that the unitary transformation is not identity and also S =/= C, then the vanishing of 
potential is due to Z = 0. Then, from (40), it is clearly found that a continuous increase of the 
coupling strength Za(S 2 — C 2 ) from zero, electron states can be pushed up in to the positive 
energy continuum. For Za(S 2 — C 2 ) <C 1 the energy formula can be expanded as: 



£ r ,„-m Z 2 a 2 (S 2 -C 2 ) 2 



m 



1 Z 2 a 2 , 1 3(5' 2 - C 2 f 
I ( 1 — - 

n 2 2n 3 m 



An 



It is seen that for Za(S 2 — C 2 ) = the modified binding energy, = E r>n — m, is 



E^ = E r , n - m = 0, 



(44) 



(45) 



where (mb) := (modified binding). Hence, we see that with increasing Z the modified binding 
energy increases also. Note that (40) is nearly real for all values of Z and n. Namely, for the 
states with n = j + 5, energy value can be calculated as: 



m 



n 2 + Z 2 a 2 



n 2 + {2CSZa) 2 



(46) 



Therefore, this results show that for all unitary transformations which apply to Dirac equation 
for gauge fixing condition and Schrodinger like requirement, the energy eigenvalues of the states 
with n = {j + \) is real for all values of Z (see the second footnote on the page 7) . 
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